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Abstract. We present new 450 and 850 fim SCUBA data of the dwarf galaxy NGC 1569. We construct the 
mid-infrared to millimeter SED of NGC 1569, using ISOCAM, ISOPHOT, IRAS, KAO, SCUBA and MAMBO 
data, and model the SED in order to explore the nature of the dust in low metallicity environments. The detailed 
modeling is performed in a self-consistent way, synthesizing the global ISRF of the galaxy using an evolutionary 
synthesis model with further constraints provided by the observed MIR ionic lines and a photoionisation model. 
Our results show that the dust properties are different in this low metallicity galaxy compared to other more 
metal rich galaxies. The results indicate a paucity of PAHs probably due to the destructive effects of the ISRF 
penetrating a clumpy environment and a size-segregation of grains where the emission is dominated by small 
grains of size ~ 3 nm, consistent with the idea of shocks having a dramatic effect on the dust properties in 
NGC 1569. A significant millimetre excess is present in the dust SED which can be explained by the presence of 
ubiquitous very cold dust (T = 5 — 7 K). This dust component accounts for 40 to 70 % of the total dust mass 
in the galaxy (1.6 — 3.4 x 10 Mg) and could be distributed in small clumps (size ~ a few pc) throughout the 
galaxy. We find a gas-to-dust mass ratio of 740 - 1600, larger than that of the Galaxy and a dust-to-metals ratio 
of 1/4 to 1/7. We generate an extinction curve for NGC 1569, consistent with the modeled dust size distribution. 
This extinction curve has relatively steep FUV rise and smaller 2175 A bump, resembling the observed extinction 
curve of some regions in the Large Magellanic Cloud. 
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1. Introduction 

The absorption of stellar radiation by dust and its subse- 
quent reemission in the infrared (IR) to the submillimeter 
(submm) is a fundamental process controlling the heating 
and cooling of the interstellar medium (ISM; e.g. Tielens 
& Hollenbach 1985). The dust IR spectral energy distri- 
bution (SED) of a galaxy is its footprint reflecting fun- 
damental physical parameters such as initial mass func- 
tion (IMF), age, stellar population and metallicity. Thus, 
knowledge of the physical characteristics of dust in galax- 
ies opens the door to understanding the star formation his- 
tory and the evolution of galaxies. However, an accurate 
interpretation of the SEDs of galaxies requires knowledge 
of the detailed macroscopic and/or microscopic dust prop- 
erties in galaxies, such as composition, emissivity, dust size 
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distribution, and spatial distribution. How these proper- 
ties are effected by variations in the metallicity of the ISM 
is the focus of our present study. 

Dwarf galaxies in our local universe are ideal labora- 
tories for studying the interplay between the ISM and 
star formation in low-metallicity environments (Hunter 
& Gallagher 1989). They are at relatively early epochs 
of their chemical evolution, possibly resembling distant 
protogalaxies in their early stages of star formation. The 
subject of dust formation in primordial environments and 
the processes controlling the subsequent evolution of dust 
compels us to study the dust properties in nearby dwarf 
galaxies, with the eventual goal of constructing templates 
that would apply to primordial galaxies. While only a 
relatively small number of metal-poor galaxies has been 
observed at mid-infrared (MIR) wavelengths using the 
Infrared Space Observatory (ISO), it is already apparent 
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that the characteristics of the MIR dust components differ 
remarkably from those of more metal-rich starburst galax- 
ies (e.g. Lutz et al. 1998; Crowther et al. 1999; Thuan et 
al. 1999; Madden 2000; Vacca et al. 2002; Madden et al. 
2003). Is it correct to simply assume Galactic dust prop- 
erties when modeling the dust in dwarf galaxies? 

Discrepancies are arising as to the quantity of dust that 
is actually present in even the lowest metallicity galaxies 
that are observed today. Since optical observations are 
more readily accessible than MIR to far infrared (FIR) 
wavelengths which require space and air borne platforms, 
the dust opacity has mostly been determined through the 
effects of obscuration, not from the dust emission prop- 
erties. These different measurements are not always con- 
sistent. Even at 1/50 solar metallicity, I Zw 18 contains a 
non-negligible amount of reddening, determined from op- 
tical observations to be equivalent to A v = 0.5 (Cannon 
et al. 2002). In the dwarf galaxy SBS 0335-052 (1/40 solar 
metallicity), A v values as high as 20 to 30 have been sug- 
gested from MIR observations (Thuan et al. 1999; Plante 
& Sauvage 2002), contrary to the low values, A v < 1, ob- 
tained by Dale et al. (2002) and Izotov et al. (1997). All 
of the details that go into modeling the dust properties in 
dwarf galaxies are critical for determining the quantity of 
dust that resides in these metal-poor galaxies. This must 
also be reconciled with current dust evolution models be- 
fore we can understand the early epoch of star formation. 

We focus on the detailed modeling of the SEDs of 
dwarf galaxies, constructing the SEDs with as much wave- 
length coverage as possible. We began with a sample of 
dwarf galaxies observed in the MIR (Madden et al. 2000; 
Madden et al. 2003). IRAS and ISOPHOT observations 
are included, along with our recently-obtained ground- 
based submm (Galliano et al. 2002) and millimetre (mm) 
observations (Galliano et al. 2003). The basis of the dust 
modeling is the model of Desert et al. (1990; hereafter 
called DBP90). To synthesise an input interstellar ra- 
diation field (ISRF), we use a stellar evolution model 
(PEGASE; Fioc & Rocca-Volmerange 1997) and addi- 
tional constraints obtained from a photoionisation model 
(CLOUDY; Ferland 1996). Here we present the complete 
modeling for NGC 1569, a nearby (D = 2.2 ± 0.6 Mpc; 
Israel 1988) dwarf galaxy with an average metallicity of 
1/4 solar (Gonzalez Delgado et al. 1997), currently in the 
aftermath of a massive burst of star formation (Israel 1988; 
Waller 1991). A dust SED model has recently been pre- 
sented by Lisenfeld et al. (2002) . Recent Chandra observa- 
tions find large inhomogeneities in the metal abundances 
with ranges from 0.1 Z to 1 Z© (Martin et al. 2002), evi- 
dence that the ISM has been affected by numerous super- 
novae explosions since the last burst. NGC 1569 contains 
two bright super-star-clusters (SSCs) (Arp & Sandage 
1985; Aloisi et al. 2001) which have blown out a large 
fraction of the gas in a dramatic display. 

The paper is organised as follows. Section [21 presents 
an overview of our new observations and the data we use 
from the literature. Section|31describes the adopted model 
and the method we use to investigate the dust properties. 



Section H presents the modeled SED of NGC 1569 and the 
consequences of the results on the dust properties. We end 
with a summary and the conclusions in Sect.0 

2. The Observations 

An observed SED was constructed, as completely as pos- 
sible, incorporating data from the literature for various 
telescopes: IRAS, KAO, ISO (ISOPHOT) and IRAM 
(MAMBO), as well as our ISOCAM data (Madden et 
al. 2003) and our new 450 and 850 [im JCMT (SCUBA) 
observations presented here. Since we are modeling the 
global SED in this paper, the differences in beam sizes are 
not of great concern here. 

2.1. SCUBA observations 

We obtained 450 and 850 [im data of NGC 1569, with 
SCUBA (Holland et al. 1999), a bolometer array on 
the James Clerk Maxwell Telescope (JCMT), during two 
observing runs in February 2000 and December 2000. 
Observations were carried out in the jiggle-mapping mode 
using a 64-point jiggle pattern with a chop throw of 150". 

We reduced these data using the SURF VI. 5-1 soft- 
ware (Jenness & Lightfoot 1998), other STARLINK utili- 
ties such as KAPPA, CONVERT, FIGARO and our own 
IDL routines. We scrutinised the data between each step 
of the processing to check for possible artifacts. We pro- 
ceeded in the following way: 

1) Atmospheric extinction: The relatively low atmospheric 
transmission is the main difficulty in obtaining good qual- 
ity ground-based observations at these wavelengths. To 
quantify the atmospheric extinction and correct our data 
for this effect, we performed skydips at 450 /im and 
850 fim every 1-2 hours. In addition, we used the mea- 
sured tcso a t 225 GHz, from the Caltech Submillimeter 
Observatory radiometer, which ranged from 0.04 to 0.1 
during our observations. We evaluate the atmospheric 
transmission at 850 /zm, rsso, using the skydip value 
adapted to the elevation of the bolometer. However, the 
450 /im skydip has been determined to be unreliable due 
to the fluctuations of the atmosphere. We, therefore, used 
relationships between tcsOj tsso and T450 to compute T450: 
t 450 = 6.52 x (t 850 - 0.049) (Archibald et al. 2000). 

2) Flatfields, despiking and bad bolometers: The non- 
uniformity of the response of each bolometer is corrected 
with standard calibration Hatfield measurements since it 
remains constant with time and does not need to be re- 
measured every night (Jenness & Lightfoot 1998). The 
spikes were removed from individual maps using a 4<r 
threshold and noise measurements were performed at the 
beginning and at the end of the shift every night. Bad 
bolometers were removed with a threshold of 3a. This last 
step can create holes in individual maps. However, since 
the array rotates on the sky, the final map contains no 
holes. 

3) Sky fluctuations: The fluctuations in the emissivity of 
the atmosphere were removed using the fact that the outer 
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ring of bolometers was observing the sky. The median of 
the outer ring, excluding the noisy bolometers, was sub- 
tracted from each map. Using cumulative growth-curves, 
we verified that no source flux was removed. Each map was 
rebinned with a pixel field of view of 6"/pixel at 850 fim 
and 3"/pixel at 450 fim which correspond to the jiggle- 
map sampling. 

4) Flux conversion: Each night, we performed three mea- 
surements of the calibration sources, Uranus, Mars or 
CRL 618, using the same observing parameters. An av- 
erage Flux Conversion Factor (FCF) was computed from 
these measurements, for each night and the individual 
maps were calibrated by the corresponding FCF. We eval- 
uated the average peak-to-aperture ratios, Rdate(A), on 
the calibrators, to correct for the loss of flux in secondary 
lobes. The full-width-half maximum (FWHM) beam size 
is 8.5" at 450 fim and 15.2" at 850 fim. 

5) Map coadding and filtering: The individual maps were 
coadded, weighted by Ti n t/<7 2 , where Tint is the integra- 
tion time and a the standard deviation of the pixels in 
the image. We applied the MR/1 multiresolution filtering 
method (Starck et al. 1998) using a 3<r threshold. 

6) Evaluation of uncertainties: To estimate the cumula- 
tive uncertainty AF U (X) on the flux Ti,(A) integrated in 
a circular aperture O (we took O = 120" for the entire 
galaxy), we quantified the individual error contributions. 
We first construct RMS maps for 850 fim and 450 /im 
observations such that each pixel A) of these maps is 
the standard deviation <trms(*7.7; A) of the corresponding 
pixels of the individual maps. The contribution of these 
fluctuations to the total error, A-F RMS (A), is computed by 
integrating the flux in the RMS map in the aperture 0: 

AF„ RMS (A) = ^£ (iJ - )ee Kms A)iWA)7V bcam ) 2 
where iVbeam is the number of pixels per beam. The error 
due to the sky emission subtraction, <r s ky, is the standard 
deviation of the distribution of the points used to com- 
pute the median: AF* k y (A) = v/V (CT sky i?datc(A) V bcam ) 2 
where Nq is the number of pixels inside the aper- 
ture 0. The error on the absolute calibration, consid- 
ering the distribution of the FCF, {C ca i(A)}, and the 
distribution of the peak-to-aperture ratios, {R ca i(A)}: 
Af™ v (A) = <7(Ccai(A))/(C C ai(A)) x F„(A), AF- tio (A) = 
<7 (-Real (A)) /(-Real (A)) x F V (X). The total uncertainty on 
the net flux is the sum of these different contributions 
(Table [TJ: 

AF V {\) = [(AF RMS (A)) 2 + (AF^ ky (A)) 2 

+ (AF„ conv (A)) 2 + (Ai^' atio (A)) 2 ] 1/2 . (1) 

7) Radio and molecular contributions: The radio contin- 
uum emission contaminating the 450 fim and 850 fim 
broadbands was evaluated from Israel & de Bruyn (1988). 
Taking their value of the spectral index, a = —0.36, we 
extrapolate the radio continuum to be 49 mJy (~ 14 %) at 
850 fim and 39 mJy (~ 3 %) at 450 fim. We estimate the 
CO (3-2) contribution in the observed 850 fim band from 
Meier et al. (2001) to be - 5 %. The observed 450 fim 
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RMS 


Sky 


Calibrator 


Ratio 


Total 


850pm 


8% 


5% 


5% 


4% 


12% 


450/im 


38% 


4% 


20% 


14% 


45% 



Table 1. The various contributions to the uncertain- 
ties on the net flux of the SCUBA maps integrated in a 
120" diameter aperture centered at a(2000) = 4 h 30 m 49 s , 
5(2000) = 64°50'55". 



and 850 fim flux values we use for our dust modeling were 
adjusted for non-dust contamination. 

Our 850 (im value (345 ± 40 mJy) is consistent with 
Lisenfeld et al. (2002) (410 ± 45 mJy) using different 
SCUBA data. Our 450 fim flux is 1320 ± 450 mJy while 
Lisenfeld et al. (2002) report a higher total flux of 1820 ± 
700 mJy but still consistent. By studying the behavior of 
both the 450 and 850 fim skydips as a function of time 
of day and sky conditions, we gave particular attention to 
calibration questions. The most accurate calibration uses 
the 850 fim skydips to adjust the 450 fim maps. Using the 
450 fim skydip to calibrate the 450 fim data is known to 
give higher fluxes (Archibald et al. 2000). This is proba- 
bly the greatest source of difference in the final calibration 
between the present paper and Lisenfeld et al. (2002). 

2.2. Infrared to millimetre data 

The ISOCAM circular variable filter (CVF) spectrum pro- 
vides the constraints in the 5 to 16 fim wavelength range. 
The details and the data treatment of the CVF obser- 
vations are presented in Madden et al. (2003). We have 
found that the slope of the MIR CVF is a critical fac- 
tor in constraining the dust model. Thus, using only the 
IRAS 12 fim band or several ISOCAM broad band obser- 
vations do not provide sufficient constraints on the model. 
We characterised the MIR dust continuum by choosing 10 
wavelength regions of the CVF which do not contain aro- 
matic bands or ionic lines (Table |5J). This was to balance 
the full observed SED which contains 8 additional points 
from other instruments. The ISOCAM images (Fig.Q) are 
2x2 raster maps and have been processed in the same way 
as the CVF spectrum (Madden et al. 2003). We decon- 
volved the individual images using a multiresolution Lucy 
algorithm from the MR/1 package (Starck et al. 1998), us- 
ing a 3<7 detection threshold. Some ratio maps are shown 
in Fig. □ 

IRAS fluxes for NGC 1569 are given by Thronson & 
Telesco (1986), Hunter et al. (1989a) and Melisse & Israel 
(1994). The values for the four IRAS broadbands given by 
Thronson & Telesco (1986) and by Hunter et al. (1989a) 
are similar, and equal to the ISOCAM CVF when inte- 
grated over the equivalent IRAS 12 fim band. The values 
reported by Melisse & Israel (1995) differ significantly and 
are flagged by the authors to be uncertain. For our mod- 
eling purposes, we use the values of Hunter et al. (1989a) 
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since error bars are also provided (Table These fluxes 
have been color-corrected. 

ISOPHOT is an imaging polarimeter operating be- 
tween 2.5 and 240 fim, which was onboard the satellite 
ISO (Lemke et al. 1996). Data were taken by the two 
cameras of PHOT-C, CU00 and C.200, at 60, 100 and 
180 fim. However after analysing the data, we realised 
that the background position was on a bright cirrus spot, 
rendering the data difficult to interpret. Instead, we use 
only the value from the ISOPHOT 170 /urn serendipity 
survey (Stickel et al. 2000) which gives a total 170 fim 
flux of (28 ± 9.8) Jy for NGC 1569 (Martin Haas, private 
communication) . 

Hunter et al. (1989b) measured a 155 fim flux of 36 Jy 
on the Kuiper Airborne Observatory (KAO). The 1.2 mm 
IRAM flux is from Lisenfeld et al. (2002). They have sub- 
tracted the contributions of the radio continuum and the 
CO(2-l) line emission. All of the fluxes used are listed in 
table |U 



Instrument 


A (fim) 


AA (fim) 


Net flux (mJy) 


ISOCAM 


8.8 


0.05 


(370 ± 220) a 




10.1 


0.96 


(480 ± 230)" 




10.8 


0.87 


(610±210) a 




11.8 


0.73 


(860 ± 220) a 




12.1 


0.69 


(980 ± 220) a 




13.2 


0.54 


(1070 ± 220) a 




13.9 


0.45 


(1270 ± 230) a 




14.6 


0.35 


(1420 ± 240) a 




15.0 


0.30 


(1470 ± 250) a 




16.0 


0.16 


(1990 ± 290) a 


IRAS 


24.3 


6.3 


(8550 ± 1100)" 




62.8 


19 


(48900 ± 7000)" 




103 


22 


(54800 ± 8000)" 


KAO 


155 


30 


(36000 ± 11000) c 


ISOPHOT 


174 


90 


(28000 ± 9800) d 


SCUBA 


443 


20 


(1280 ± 450) e 




863 


70 


(280 ± 60) e 


IRAM 


1200 


200 


(190 ±60) ; 



Table 2. Net fluxes for the dust emission by the entire 
galaxy used for the modeling. The values take account of 
all the corrections mentioned in the text. ° Madden et al. 
(2003), b Hunter et al. (1989a), c Hunter et al. (1989b), 
d Stickel et al. (2000) & full catalog in prep., e this paper, 
f Lisenfeld et al. (2002). 



2.3. UV and optical data 

To constrain the input stellar radiation field (Sect. l3~3*)l we 
used optical and UV data measured for the total galaxy 
(Israel 1988; De Vaucouleurs et al. 1991; Tabled)). 



Wavelength 

(A) 


Flux density 

(10- 14 erg s- 1 cm 1 A" 1 ) 


1500 


(152 ± 48) a 


1800 


(208 ± 25) 


2200 


(146±25) a 


2500 


(100 ± 18) a 


3650 
4400 
5500 


(108 ± 10)" 
(101 ±9)" 
(70 ± 6)" 



Table 3. Fluxes for the entire galaxy used to model the 
global ISRF. The values are corrected for Galactic extinc- 
tion. a Israel (1988), 6 De Vaucouleurs et al. (1991). 



3. Self-consistent modeling of the global SED 

We have modeled the IR to millimetre dust emission from 
NGC 1569 in a self-consistent way using a wide variety 
of observational constraints. In this section we describe 
the method used. The results are presented in Sect. 0J 
The general algorithm used to obtain the best solution is 
outlined in Sec. 13.41 

3.1. Physical parameters of the galaxy 

In the adopted model, we assume that the dust is located 
in a thin shell at an effective distance (R e ff) from the cen- 
tral heating sources, the stars (Fig. |HJ • The extinction is 
treated in terms of a slab model which neglects scattering. 




Fig. 3. Schematic diagram of the model geometry. 
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Fig. 1. ISOCAM and SCUBA images of NGC 1569. Top-left: ISOCAM-LW6 (A = 7.75 fitn, beam=7.6"), this band 
traces PAH features. Top-right: ISOCAM- LW3 (A = 15 fim, beam=9.9"), this band traces the hot dust continuum. 
Bottom-left: SCUBA (A = 850 fim, beam=15.2"). Bottom-right: SCUBA (A = 450 fim, beam=8.5"). These two 
submm bands trace the cold dust continuum. The field of view is the same for the 4 images and the color table has 
the same dynamic range. The ISOCAM images have been deconvolved. 



"Kit: lJtJ£? ''liin^o ^tnj.-'d^C miil (ji.mUi.Lt l:j <in--i Mif.' I'jiii: ■. Lljj l^l^l - : 1. ^tnj.-'d^C urn <.i:i'A (ji.mUi.Lt T'.vc ,uiu'} .;-.!■>. r>+ 




si-" as? oi- 1 aj- u cc^ :o v 



Fig. 2. Ratio of hot dust/cold dust. Left: ratio of the ISOCAM-LW3 image (contours, A = 15 jim) to the SCUBA image 
(A = 850 fim). Right: ratio of the ISOCAM-LW6 image (contours, A = 7.75 fim) to the SCUBA image (A = 850 {im). 
The ISOCAM images have been degraded to the resolution of the SCUBA (850 fim) image (15"). 
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In order to have a consistent set of physical parameters 
(average radius, density, column density, etc.), we estimate 
the radius R oqu iv of the equivalent spherical galaxy from 
our 850 fim image and deduce the value of the densities 
from this radius and from the hydrogen masses given in 
the literature. To determine R oqu iv, we fit an ellipse to the 
galaxy at 850 fxm using the geometrical mean of the axes 
as the average radius, which gives an equivalent radius, 
Rcquiv, of 0.65±0.18 kpc (the uncertainty, AR cq uiv, is due 
to the uncertainty in the distance), consistent with the 
optical size: 1.85 x 0.95 kpc (Israel 1988) which gives a 
Rcquiv = 0.66 kpc. The corresponding effective hydrogen 
density, n ff, and column density, N e fj, are 6.7 cm~ 3 and 
1.8 x 10 22 cm" 2 from the M(HI) and M(H 2 ) given by Israel 
(1997) (M(HI) = 1.4 x 10 8 M , M(H 2 ) = 0.5 x 10 8 M ). 

The mean distance from the dust to the stars is taken 
to be R e ff = 3/4 x R oqu i v = 0.50 kpc, corresponding to a 
homogeneous spatial distribution of the dust. We ran the 
model for this value of R c ff and also for two extreme cases: 
R™ ff ax = Rcquiv + AR equiv = 0.84 kpc and R™ n = 1/2 x 
Rcquiv — AR cqu i V = 0.24 kpc. These extreme values permit 
us to investigate a wide variety of spatial distributions and 
also to take into account uncertainties on the distance 
determination. 

3.2. Dust modeling 

To compute the dust emission spectrum we use DBP90 
which includes three dust components and provides a co- 
herent interpretation of both the interstellar extinction 
and the infrared emission. This model has been used to 
explain the dust properties in the Galaxy. The three com- 
ponents of the DBP90 model are the following: 

1. Polycyclic Aromatic Hydrocarbons (PAHs); 2- 
dimensional molecules responsible for the MIR 
emission features and the FUV non-linear rise in 
the extinction curve. In the standard DBP90 model 
applied to the Galaxy, PAHs have a size range of 
0.4 nm < a < 1.2 nm. 

2. Very Small Grains (VSGs); 3-dimensional carbona- 
ceous grains responsible for the MIR continuum emis- 
sion and the 2175 A extinction bump. In the standard 
model, their radii are 1.2 nm < a < 15 nm. 

3. Big Grains (BGs); 3-dimensional silicates, coated or 
mixed with carbonaceous material, responsible for the 
FIR emission and for the NIR and visible rise of the 
extinction curve. In the standard model, their radii are 
15 nm < a < 110 nm. 

Due to their small sizes, and thus low heat capacity, PAHs 
and VSGs are transiently heated by single photon absorp- 
tion and are not in thermal equilibrium with the radiation 
field. BGs usually reach thermal equilibrium for the grain 
sizes determined for the Galaxy. However the dust model 
does include thermal fluctuations for the BGs if their sizes 
warrant it. The model also computes the absorption cross 
section of the grains corresponding to the physical param- 
eters that fit the emission. Finally, the extinction curve 



is generated using the absorption cross sections generated 
by the model. 

We use the DBP90 model to compute the MIR to mil- 
limetre SED. The data which constrain the SED are sum- 
marised in TableEl To obtain the best fit to these data, we 
have varied the physical parameters of the model, using a 
Levenberg-Marquardt method (Numerical Recipes, Press 
et al. 1996) to minimize the \ 2 - The parameters which can 
vary are: 

1. the mass abundance of each species, Y = m/mn (m 
is the dust mass in the beam and wh is the hydrogen 
mass in the beam); 

2. the index of the power-law of the size distribution of 
each species, a, where n(a) cx a~ Q (a is the grain radius 
and n(a) the number density of grains between a and 
a + da); 

3. the minimum and maximum sizes of each species, a m i n 

and dmax ) 

4. the absorption coefficient of the grains, Q a hs, defined as 
Cabs = Qabs^a 2 (cr a bs is the absorption cross-section) 
for a spherical grain of radius a. 

Moreover, the widths of the individual broadband obser- 
vations are taken into account, except for the ISOCAM 
CVF spectrum since the spectral resolution is sufficient 
to neglect this effect (R > 35). For the data from the 
bolometers (SCUBA and MAMBO), we integrated the 
modeled SED into the observational broadbands. For the 
data from photo-multipliers (KAO and ISOPHOT), we 
made the corresponding color-correction when integrating 
into the broadbands. For the IRAS data, this correction 
was already made. This correction has the largest effect 
for the 170 ISOPHOT broadband (Fig.©. 

3.3. Modeling of the interstellar radiation field 

The DBP90 model requires an input ISRF which heats the 
dust. To be consistent, we have modeled the global ISRF 
for NGC 1569 using the UV and optical data (Table EJ 
coupled with the evolutionary synthesis model PEGASE 
2.0 (Fioc & Rocca-Volmerange 1997), taking into ac- 
count the constraints from photonisation processes using 
CLOUDY. 

3.3.1. Stellar evolutionary synthesis 

PEGASE computes the stellar spectral energy distribu- 
tion taking into account the metallicity evolution. We do 
not incorporate the dust extinction effects and the nebular 
emission into PEGASE. The Q a b s values are computed by 
DBP90, the extinction is calculated with a simple screen 
model and the nebular emission is taken into account us- 
ing the photoionisation model, CLOUDY fSect. l3~3~2"|) . 

We fit the UV-to-optical data (Table |3J) using a com- 
bination of two single instantaneous bursts with an initial 
metallicity of Z = Z Q /4.27 (Gonzalez Delgado et al. 1997). 
The reddening is calculated using the dust cross-sections 



F. Galliano et al.: The Dust Spectral Energy Distribution of NGC 1569 



7 



computed by the dust model. However, the optical depth 
deduced from this extinction curve corresponds to the case 
where all the dust is located in front of the stars. To obtain 
the effective optical depth, r® ff , the shape of the extinction 
curve is scaled in order to satisfy the energy conservation: 



F* — FrjV-opt + FiR_ r 



(2) 



where F* is the flux density as it would be observed if 
there were no dust in the galaxy, only stars, and Frjv-opt 
and FiR_ min are the flux densities from reddened stars 
and dust respectively, as they would be observed: 



^UV-opt _ 
FrjV-opt = 
FlR- mm = 



F* exp(-r, 



eff v 



F*exp(-rf)di/ 
F* [l-exp(-rf)] Av. 

-OO 



(3) 



We minimise the \ 2 to find the best age combination. 
The fit gives several age combinations with similar x 2 ■ To 
remove this degeneracy, we use a further constraint, com- 
paring the MIR ionic line ratios from Madden et al. (2003) 
and the theoretical results computed with the photoioni- 
sation model CLOUDY (version 90.04, Ferland 1996). 

3.3.2. Photoionisation 

CLOUDY predicts the spectra of astrophysical plasma in 
different environments. It uses a recent atomic database 
and takes into account the geometry and the gas proper- 
ties. We use CLOUDY through the MICE 1 IDL interface. 
A thick shell geometry is adopted while that for the dust 
emission is a thin shell. 

We have computed a grid of solutions from CLOUDY 
containing a range of values for the following parameters: 

1. the shape of the input ISRFs are those which give the 
best x 2 values as the output from PEGASE; 

2. the total luminosity of the galaxy was computed by 
integrating the total UV-optical and IR-mm SED: 
L tot = 1.9 x 10 9 L Q ; 

3. the inner radius varies from Rj n — 5 pc to Rj n = 90 pc, 
which is the radius of the hole in HI centered on the 
SSC A (Israel & Van Driel 1990); 

4. the average hydrogen number density is n(H) e fi 
(6.7 cm" 3 ; Sect. ED; 

5. the outer radius was taken to be R oq uiv CSect. l3~TJl . 

We set the elemental abundances to X = Yjj = X®, Y = 
Y Ho = Y Q and Z = 1 - X - Y = Zq/4.27 according to the 
value of the metallicity given by Gonzalez Delgado et al. 
(1997). Fitting the theoretical ionic line ratios to the MIR 
data ([NeIII]/[NeII] = 9.6±1.4; [SIV]/[NeIII] = 0.49±0.08; 
Madden et al. 2003) clearly removes the degeneracy of the 
ISRF. 



1 MICE is supported by the SWS and the ISO Spectrometer 
Data Center at MPE through DLR (DARA) under grants 50 
QI 86108 and 50 QI 94023. 



3.4. Iterative process 

The modeling of the dust emissivity generates an extinc- 
tion curve which is used to deredden the UV-to-optical 
data in order to compute the ISRF. The ISRF is used by 
the dust emissivity model to heat the dust. Therefore, to 
be consistent, we compute this sequence through an iter- 
ative process until we reach a stable, self-consistent solu- 
tion. This solution is obtained when the extinction curve 
used to deredden the data equals that deduced from the 
dust emission. The general scheme of this iterative process 
is shown in Fig. 0] 

UV-opt. Obs. 



PEGASE 






1 


I 



Several 
ISRF 



i 



Extinction 



IR-mm Obs 




IR ionic lines 



IR-mm SED 



Fig. 4. Algorithm used to compute the SED in a self- 
consistent way. "UV-opt. Obs." are the data contained in 
Tabled "IR-mm Obs." are the data contained in Tabled 
"IR ionic lines" are from Sect. 13.3.21 "IR-mm SED" is 
the modeled dust SED and "Extinction" is the extinction 
curve computed by DBP90 after converging to a solution. 
We iterate on this process to obtain the necessary consis- 
tency between emission and extinction. 



4. Results and discussion 

We have applied the general model described in Sect. El 
to our observed IR-to- millimetre SED for NGC 1569 
(Table [2J. 



4.1. The dust emissivity 

In this subsection we describe the details of the solution 
of the dust modeling and discuss the constraints and the 
consequences of the model results. 



4.1.1. Millimetre excess 

In areas where we have no strong constraints, some param- 
eters of the general model described in Sect. 13.21 are kept 
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fixed, while others are left as free parameters to obtain 
the best fit. 

1. PAHs: this component is too weak in NGC 1569 to be 
constrained with our set of data. An upper limit on 
the mass abundance of the PAHs, Ypah, is fixed and 
we used the Galactic values for a, a m j n and a max . 

2. VSGs: we varied the mass abundance, lySGi an d two 
parameters governing the size distribution, a and a max . 
However, we cannot constrain a m ;„ since the \ 2 P ro ~ 
cess tends to decrease it to very small sizes which have 
no physical meaning, so it was kept constant at the 
Galactic value. 

3. BGs: we varied the mass abundance, YbG; & n d two 
parameters governing the size distribution, a and a m i n . 
Since the results are not sensitive to a max , we adopt 
the Galactic value. 

We are able to fit the SED quite well except in the mil- 
limetre wavelength range, where there is an excess that 
can not be explained by the three standard DBP90 com- 
ponents (see Fig. EJ . 

Although very cold dust is not the only likely expla- 
nation of this excess, we have added a fourth component 
modeled with a modified black-body to describe it. We 
call this component Very Cold Grains (VCGs). We dis- 
cuss the other possible origins and the consequences of 
this excess in Sect. 14.31 The monochromatic flux radiated 
by the VCG component in a sphere of radius R c ff is given 
by: 

VCG = 3Y VCG M(H) /Qo\ fXoV 2hc/A 3 

ipR 2 eS \ a J \ A J exp(hc/AkT) - 1 

(4) 

where M(H) is the total mass of hydrogen in the galaxy, 
Ao, Qo and p are the reference wavelength, the absorp- 
tion coefficient at Ao and the mass density of the grains. 
This VCG radiated flux depends on three variables: the 
temperature, T, the emissivity index, (3, and the mass 
abundance of the grains, Ivcg (we adopted typical val- 
ues: Q /a = 40 cm -1 , A = 250 pm, p = Psiiicatc = 
3.3 gcm~ 3 if fi = 2 and p = p gra phito = 2.2 gcm~ 3 if 
/3 = 1). As we have only three data points to constrain the 
submillimetre-millimetre wavelengths regime of the SED 
(450 pm, 850 pm and 1.2 mm), we consider a way to fix the 
value of one parameter to avoid the degeneracy. Since the 
best solution requires a relatively flat slope, we fix (3 — 1 
for the x 2 solution. Setting (3 = 2 gives us an unphysical 
solution with a maximum temperature of 3 K. Our emis- 
sivity index of 1 is lower than that found by Dunne & 
Eales (2001) for a sample of IR-bright galaxies ((3 ~ 2). 
However, the assumption of a lower (3 is reasonable in light 
of the fact that Agladze et al. (1996) measure the absorp- 
tion coefficient of cosmic dust analog grains and find that 
over the temperature range of 1.2-10 K, (3 decreases with 
temperature. Thus, the mass abundance, YycGj & n d tem- 
perature, T, remain the only free parameters governing 
the VCG properties. 



4.1.2. What constraints can we put on the dust 
properties ? 

Our best fit of the observed SED is shown in Fig. |SJ The 
values of the parameters corresponding to the best \ 2 fit 
for the 4 component model arc summarized in Tabic 
This corresponds to \ 2 — 1-5 with 9 degrees of freedom. 
In order to quantify the reliability of the grain proper- 
ties deduced from our model, we attempt to estimate, as 
conservatively as possible, the errors on these parameters. 

We investigated the two possible sources of uncertain- 
ties: the sensitivity of the model itself to certain parame- 
ters, and to the value of the average radius (R c ff ) which in- 
fluences the radiation field intensity. First, we estimate the 
errors on the goodness of the fit, by quantifying the depar- 
ture of the model from the observations. We independently 
vary each parameter to fit the extreme cases allowed by 
the error bars. For example, the VSGs are constrained 
mainly by the observations from 8 to 16 pm. Thus, we 
made 4 extreme fits of the VSGs: one through the upper 
edges of the error bars, one through the lower edges, one 
through the upper edge of the error bar at 8 pm and the 
lower edge at 16 pm, and one through the lower edge at 
8 pm and the upper edge at 16 pm. This gives us an idea 
of the sensitivity of the model to the physical properties. 
Second, the model was solved for the two extreme values 
of the radius (R^ ax and R^ n defined in Sect. ED . These 
two extreme cases give us an idea of the influence of the 
geometry (and also of the ISRF) on the dust properties. 
Figure [7| demonstrates the effect of the assumed radius on 
the grain size distributions. 



10 c r R = R m „ 

: R = R m „„ I 

: R = R„„„ i 




1 10 100 

a (nm) 



Fig. 7. Grain size distribution for the two extreme radius 
cases (R£fj x and RJ,ff n ) and for the best estimated value, 
R c ff (see section l3~TJl . It gives the conservative range within 
which the size distribution can be constrained. 



Table[5]contains the ranges of the physical parameters 
for the dust components given the sources of uncertainty 
from the fitting and the geometry. The PAH and VSG dust 
parameters are not very sensitive to the assumed geometry 
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Fig. 5. NGC 1569 observations and modeled SED. The data (Table 0) are indicated by crosses: vertical bars are the 
errors on the flux values and the horizontal bars indicate the widths of the broadbands. The lines are the dust model 
and its different components (see Sect. 14. lTTfl . Diamonds are the model integrated over the observational broadbands 
and color-corrected. Thus, departures from the model lie where the diamonds deviate from the crosses. The power is 
expressed in WH" 1 which is vL u divided by the number of H atoms (h u is the monochromatic luminosity). 
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Table 4. The values are those for the best fit. For comparison, we give the corresponding values for the Galaxy from 
Desert et al. (1990). 



but are relatively sensitive to the fitting. For these cases, 
Umax can change, at most, by a factor of 3. The best esti- 
mate R e ff shows a flat size distribution for the VSGs, and 
a small slope for the 2 extreme geometry cases (Fig. 0). 
The slope of the BG size distribution, however, is very 
sensitive to the assumed geometry If R e ff is increased to 
R^ 35 , the size distribution shows an approximately con- 
stant size of ~ 3 nm (shown as the near vertical dashed 
line in Fig. EJ), with a deficit of larger grains. When the 
dust is located farther from the stars, the size of the grains 



must be smaller to reach the same temperature, which is 
constrained by the observations. 

We treat the sensitivity of the parameters of the VCGs 
separately. Indeed, the temperature of the VCGs is a very 
important parameter, which has major consequences for 
the derived dust mass, and is only constrained by two data 
points. In this case, the temperature is strongly dependent 
on the emissivity index, however we excluded (3 = 2 in 
Sect. l4.14l to avoid unphysical solutions. The temperature 
does not depend on the radiation field since we added 
this fourth component as a modified blackbody. The best 
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v (Hz) 

Fig. 6. Global synthesized SED for NGC 1569. The solid line is the total energy of the galaxy: the sum of the dust 
(dashed line) plus the emerging stellar radiation (dashed-dotted line). For comparison, the dotted line shows the 
non-extincted stellar radiation. 







Lower limit 


Upper limit 


PAH 


Y 





1.0 x 10" b 


VSG 


a 


2.6 


5.2 




a+ 


35 A 


120 A 




Y 


1.3 x 1(T 5 


2.3 x 10" 5 


BG 


a 


5.8 


35 




a_ 


21 A 


29 A 




Y 


3.5 x 1(T 4 


4.7 x 10~ 4 



Table 5. Range of reliability of the physical parameters 
of the grains (from errors in the fitting and variations in 
the assumed geometry). 



the lack of PAHs. Moreover, the VSG component in the 
Galaxy and in NGC 1569 emit roughly similar energies, as 
do the BG components in both galaxies. However, this en- 
ergy is emitted at shorter wavelengths in NGC 1569 than 
in the Galaxy since these grains are in a more energetic 
environment. On the long- wavelength side of the SED, the 
energy emitted from the Galaxy is ~ 10 times higher than 
that of NGC 1569. Even though the excess in the submm 
does not contribute much to the global energy budget of 
NGC 1569, most of the dust mass resides in this excess 
while in the Galaxy most of the dust mass is in the BG 
component. 



fits within 1 a variation give a range of temperatures: 
5 K < T < 7 K. The abundances corresponding to this 
temperature range are 0.4 x 10~ 3 < Yvcg iS 1-3 x 10~ 3 . 

4.1.3. The dust properties - what are the 
consequences? 

The dust properties deduced from this study differ from 
those of normal, more metal rich galaxies like the Milky 
Way. If we compare the shape of the dust SED of 
NGC 1569 to the SED of the Galaxy (Fig.®, we notice 



The dust size distribution and extinction curve 

These differences noted above are reflected in the size dis- 
tribution of the grains. Figure compares the mass spec- 
trum of the best solution of the 3 DBP90 dust components 
in NGC 1569 with that of the Galaxy. One significant dif- 
ference is that the BGs in NGC 1569, are smaller overall 
compared to those of the Galaxy. The dust mass is con- 
centrated in grains of smaller sizes in NGC 1569, ~ 3 nm, 
and the dust mass in large grains is almost negligible. 
In this case, what we called Very Small Grains and Big 
Grains have similar sizes but these names are maintained 




Fig. 8. Comparison of the dust SED for NGC 1569 (black) with the Galaxy (grey, from DBP90). In addition to the 
total dust SED (solid lines), the individual components of the DBP90 dust model also shown for both the Galaxy and 
NGC 1569. The power is expressed in WH" 1 which is vL v divided by the number of H atoms (L„ is the monochromatic 
luminosity) . 



to be consistent with the original DBP90 model. However, 
we should call these two components carbonaceous grains 
and silicate grains respectively. 

The size distribution of grains in NGC 1569 is con- 
sistent with that of grains which have been affected by 
shocks. FigurelTUlillustrates the effect of shocks on the size 
distribution of the carbon and silicate grains in NGC 1569 
(Jones et al. 1996). Both distributions have the same shape 
even though quantitatively, the cut-off does not occur at 
the same grain size. The disagreement of the shock model 
with the precise shape of the NGC 1569 dust size distri- 
bution reflects our ignorance of the initial size distribu- 
tion in NGC 1569, but the effect is qualitatively the same. 
The dust mass is transferred from large to small grains 
via shattering in grain-grain collisions in shock waves. 
NGC 1569 has recently experienced a large number of su- 
pernovae (Israel & de Bruyn 1988; Waller 1991; Greve et 
al. 2002), thus supernova-generated shock waves should 
be common in the ISM of this galaxy. 

The presence of smaller grain size distributions, that 
dominate the global IR emission in dwarf galaxies, is op- 
posite to the deduced dust properties in AGNs (e.g. Laor 
& Draine 1993; Maiolino et al. 2001). In the case of AGNs, 
dust destruction effects due to sublimation or thermal 



sputtering can explain the absence of small grains in the 
hostile AGN environment, leaving only the larger, more 
robust grains. Although similar effects could be present in 
the hard radiation field of NGC 1569, the effects of shocks 
in redistributing the dust size distribution appears to be 
dominant. However, most of the dust mass is in the VCG 
component and we cannot yet deduce the size distribution 
of these grains, until such time as the submm/mm wave- 
length range of the SED is observed in more detail. We 
give a rough estimate of the average size of the VCGs in 
Sect. EH 

Due to the small grain sizes in the PAH, VSG and BG 
components, most of the grains are stochastically heated 
and are not in thermal equilibrium with the radiation field. 
Table El contains the temperature range of each compo- 
nent for the maximum and minimum sizes compared to 
Galactic values. Stochastic heating is a very important 
process that must be taken into consideration in dust mod- 
els applied to dwarf galaxies. Using models which assume 
modified black bodies for dust emission even up to FIR 
wavelengths in dwarf galaxies, when the grains can be 
thermally fluctuating, can result in an underestimate of 
the dust mass. We have compared the grain radii at the 
transition between thermal equilibrium and stochastically 
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Table 6. Temperatures of the grains (other than VCGs). 
We give the minimum (T m i n ) and maximum (T max ) tem- 
peratures for the minimum and maximum sizes, com- 
ponent by component. All the grains are stochastically 
heated except the largest BGs where T m i n = T max = 

Tequilibrium • 



Fig. 9. Size distribution of the grains in NGC 1569 and in 
the Galaxy. The three thin lines are the mass spectrum of 
the three dust components of the Galaxy. The bold lines 
are the mass spectrum of the three dust components of 
NGC 1569. a is the grain radius and m(a) is the dust 
mass of the component between a and a + da. 
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Fig. 10. The effects of shocks on the dust size distribution. 
The solid lines are the mass spectrum of the VSGs which 
are carbonaceous grains and of the BGs which are silicates, 
in NGC 1569. The dotted lines show the mass spectra of 
carbon grains before (straight line) and after (curved line) 
a shock wave (T4 noc k = 100 kms" 1 ). The dashed lines 
represents the effect of the 100 kms" 1 shock on silicate 
grains (before shock: straight line and after shock: curved 
line). 



radiating dust (a — at) in NGC 1569 and in the Galaxy. 
First, we compute the average energy of a photon, e^(a), 

poo 

\ u„Q u (a)dv 
&r(o) = -joS , (5) 

/ (u v /hv)Q v {a) dv 
Jo 



where ti„ = 4arl v /c, Q^(a) is the absorption efficiency of 
a grain of radius a and c is the speed of light. Second, 
we calculate the temperature of the dust, Td, following a 
single absorption of a photon of energy e 7 (a) : 

47T0 3 f Td 

£r(a) = -y- J C(T)dT, (6) 

using the expression for the specific heat capacity, C(T), 
for silicate grains given by Draine & Anderson (1985). The 
transition occurs for radii where the cooling time of the 
dust, r coo i, is equal to the photon absorption time, r a bs, 
which can be estimated as follows: 

f°° kI v 
Tabs - 4?ra 2 / — ^-Q u (a) dv 
Jo e 7 ( a ) 

(7) 

1 cool — 47ra / — _ Quia) dv 
I Jo £7(0) 

7j/ is the intensity from Fig. 1151 r coo i and r a b s are plotted 
in Fig.HTIfor NGC 1569 and the Galaxy. We find that the 
transition occurs at at — 5 nm in NGC 1569 and a t ~ 9 nm 
in the Galaxy. In other words, grains larger than a*, for 
a given ISRF, will be in thermal equilibrium with the ra- 
diation field and smaller than a tl will be stochastically 
heated. Alternatively, we can show that the order of mag- 
nitude of at in the Galaxy is roughly twice at in NGC 1569. 
If we assume that €y is the same in both cases and that we 
are in the regime where Q y {a) /a is independent of a, then 
we see that r a b s oc a 3 LisRF (£isrf being the integrated 
luminosity of the ISRF) and that r coo i depends only on a. 
From the radiation fields in Fig. 1151 we know that Lisrf is 
roughly 10 times larger in NGC 1569 than in the Galaxy. 
Consequently, we expect that the transition radius in the 
Galaxy would be 10 1 / 3 ~ 3 times the transition radius in 
NGC 1569, which is the correct order of magnitude. This 
result means that grains are stochastically heated at larger 
sizes in the Galaxy than in NGC 1569. However, the size 
distributions are different. Thus, in NGC 1569, most of 
the grains are stochastically heated since the dust mass 
(PAH, VSG and BG) is concentrated in small sizes and 
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suits in a relatively steeper extinction curve. The bump at 
2175 A is weaker in NGC 1569 than in the Galaxy, due to 
the dominance of the BGs over the VSG component. This 
peculiarity is also observed in the Large Magellanic Cloud 
(Fitzpatrick 1985) extinction curve which has a metallicity 
similar to NGC 1569. Figure IT31 compares the synthesized 
extinction curve for NGC 1569 to the observed extinction 
curves for the LMC (Koornneef & Code 1981; Nandy et 
al. 1981), the SMC (Prevot et al. 1984) and the Milky 
Way (Seaton 1979). The similarity of the modeled ex- 
tinction curve with that of the LMC is striking. However, 
the intensity of the bump is, of course, model dependent. 
Nevertheless, this bump should be weaker compared to the 
continuum since the slope of this continuum is steeper. 



Fig. 11. Cooling and photon absorption times computed 
using Eq. The solid lines are for NGC 1569 and the 
dashed lines for the Galaxy. 



in the Galaxy, most of the emission originates in grains 
reaching thermal equilibrium since the dust mass is con- 
centrated in large grains. 

As a result of the difference in the grain size distribu- 
tion, the extinction curve of NGC 1569 is noticeably dif- 
ferent from that of the Galaxy (Fig. 1121) . This extinction 
curve consists only of the three standard DBP90 compo- 
nents and does not take into account the extinction due 
to the VCGs since we do not know the size distribution of 
this dust component. However, the energy emitted by the 
VCGs is very weak (~ 0.02 % of the total energy radiated 
by the dust from MIR to mm wavelengths) thus the ex- 
tinction due to this component should be very weak also. 
The magnitude of the extinction is lower in NGC 1569 
than the Galaxy at all wavelengths. Our iterative process 
gives A v — 0.45 ± 0.05. This Ay corresponds to the case 
where all the dust mass of the three DBP90 model dust 
components is located in front of the stars. The effective 
Ay deduced from the energy balance, is A^ = 0.2, less 
than the value Ay = 0.65 ± 0.04, deduced by Devost et 
al. (1997) from the Ha/H/3 line ratio toward the main 
body of the galaxy. However, they assumed Galactic ex- 
tinction properties to derive this value. Figure IT2"1 clearly 
shows that if we normalise the extinction curves of the 
Galaxy and NGC 1569 at the V wavelength, the UV slopes 
are very different indicating higher energy absorption in 
NGC 1569. Thus, to obtain the same energy absorption, 
the Ay adopted with the Galactic properties should be 
higher, which is consistent with Devost et al. (1997). In a 
general way, using Galactic extinction properties instead 
of the synthesized extinction curve of NGC 1569, when 
both are normalised by Ay, gives an erroneously lower 
t\jv by a factor as low as ~ 0.7. 

We notice in Fig. ^] that the lack of PAHs induces a 
UV rise in the extinction curve that is quite linear, due 
to the dominant effect of the BGs in the FUV rise and 
that the small size of the more abundant grains (BGs) re- 
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Fig. 13. Comparison of various extinction curves. The 
Galactic curve is from Seaton (1979), the LMC's is the 
average of Koornneef & Code (1981) and of Nandy et al. 
(1981) and the SMC's is from Prevot et al. (1984). 



The dust mass 

We determine the dust mass of NGC 1569 to be (1.6 — 
3.4) x 10 5 Mg, 15 to 20 times more than that found previ- 
ously by Lisenfeld & Ferrara (1998), who only had IRAS 
data available, and ~ 4 times more than that found by 
Lisenfeld et al. (2002) who use the DBP90 dust model, but 
assume Galactic dust grain properties and the Galactic 
ISRF. Our results show that the BGs and VCGs dominate 
the dust mass, with very little contribution from the VSGs 
and PAHs, 40 to 70 % of the total dust mass resides in 
the VCG component. The total gas-to-dust mass ratio, TZ, 
inferred for NGC 1569, assuming a total mass of hydrogen 
M H = 1.9 x 10 8 M Q (Israel 1997) and M Ho = 0.25M gas , 
is between K = 740 and K = 1600 (Table 0). If we as- 
sume only that TZ cx Z, we should find, for NGC 1569, 
TZ = 770 which is in agreement with our best % 2 value 
of 1Z = 740. However if we use the law from Lisenfeld & 
Ferrara (1998), 1/1Z oc Z a52±a25 , we find 1Z = 2900 and 




Fig. 12. Extinction curve for NGC 1569 (black) compared to the Galaxy (grey, from DBP90). The Galactic extinction 
curve has been scaled down for comparison to that of NGC 1569. The solid lines are the total extinction curve modeled 
with DBP90. The contributions to the extinction curves from the individual dust components are also shown in dashed 
lines (PAHs), dotted lines (VSGs) and dashed-dotted lines (BGs). The opacity is expressed as the cross section per H 
atom. 
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Table 7. Dust masses, gas-to-dust and dust-to-metal mass ratios in NGC 1569 deduced from the parameters in 
Table 0] The entries corresponding to non-relevant quantities are blank for more clarity. 



a lower limit of 1Z = 1200 which would be correct if we 
did not have very cold dust. It is not surprising since the 
model of Lisenfeld & Ferrara (1998) does not take into 
account submillimeter data so they did not include very 
cold dust. We also deduce the dust-to-metal mass ratio 
V = 1/(1 + KZ) which is V ~ 1/3 in the Galaxy and 
1/4 < V < 1/7 in NGC 1569. The smaller V value found 
in NGC 1569 could reflect the fact that shocks can erode 
and destroy the dust, transferring more metals into the 
gas phase. 



The spatial distribution of the dust from the observa- 
tions 

Our multi-wavelength maps can be used to study the spa- 
tial distribution of the various dust components. At first 
glance, the morphology of the submm/mm maps look 
similar to that of the ISOCAM-LW3 15 and LW6 
7.7 ^m (see Fig. QJ. Lisenfeld et al. (2002) conclude that 
in NGC 1569 both the MIR and submillimetre/millimetre 
emission are due to the same dust component. Figure [2 
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shows the spatial distribution of the ratio of LW6 (which 
traces the PAHs) to SCUBA-850 fim (which traces the 
cold dust, both BG and VCG) and the ratio of LW3 (which 
traces the VSGs) to SCUBA-850 fim (cold dust). These 
ratio maps do not show a flat distribution. They demon- 
strate that the distribution of the hot dust emission from 
the VSGs is definitely more concentrated toward the re- 
gion of the two most prominent Ha peaks whereas the cold 
dust emission and the PAHs are more extended. The PAH 
emission is more extended than the VSG's but less so than 
the cold dust emission. This is consistent with the simple 
view of an HII region and a surrounding photodissocia- 
tion region/molecular cloud. Hot dust emitting at 15 /im 
is seen to peak inside HII regions in the Galaxy, while 
the PAHs peak at the photodissociated edges of molecu- 
lar clouds (e.g. Cesarsky et al. 1996; Abergel et al. 2000; 
Klein et al. 1999). As the colder dust will be excited by 
similar sources, the morphology of the sub mm maps will 
appear to resemble the star formation tracers at these 
spatial scales. In addition, the lower metallicity ISM of 
NGC 1569, tends to be clumpier, due to the lower dust 
abundance and, hence, larger mean-free path of photons 
and, subsequently, larger photodissociation effects. This 
clumpy effect means that the pervasive non-ionising radi- 
ation will be effective at exciting the surrounding colder 
dust. We find that, contrary to commonly-made assump- 
tions, the bulk of cold dust in dwarf galaxies is not nec- 
essarily concentrated toward the outer regions, but is dis- 
tributed between the star formation sites. 



4.2. The synthesized ISRF 

The best synthesized global ISRF for NGC 1569 (Fig. EJ) 
is a combination of two single-bursts with ages 4 and 
100 Myr. The young component represents the recent 
starburst which has produced the two super-star-clusters 
(SSCs). These SSCs have populations of ages 4 — 7 Myr 
(Origlia et al. 2001), consistent with our younger popula- 
tion. The older component we use represents the underly- 
ing population present before the more recent starburst. 
Although other older populations might be present, they 
do not dominate the global ISRF for this galaxy. 

Compared to the Galaxy, the global ISRF modeled for 
NGC 1569 is much harder in the UV range due to the 
contribution from the very young population (figure El- 
At 0.1 /im, the ISRF from NGC 1569 is higher by more 
than an order of magnitude compared to the Galaxy. The 
ISRF of the Galaxy and other spiral galaxies is dominated 
by the cooler disk material while in NGC 1569, and other 
starbursting dwarfs, it is the star formation activity that 
dominates the global ISRF. This is also evident from the 
global MIR characteristics (Madden et al. 2003) and is also 
the case for other dwarf galaxies (Galliano et al. 2003). 
Therefore, when modeling the dust emission from dwarf 
galaxies, it is not valid to assume the Galactic ISRF. 




Fig. 14. Synthesized ISRF for NGC 1569 computed with 
PEGASE and CLOUDY. The Points with error bars are 
the observational data from Table 03 the solid black line 
is the global non-extincted ISRF, the dashed line is the 
global extincted ISRF and the grey line is the young 
single-burst component. The extinction curve used is the 
output from the dust model DBP90. 




0.1 1.0 



Fig. 15. Synthesized ISRF for NGC 1569 (solid line) com- 
pared to the Galaxy (dashed line). 



4.3. An explanation for the millimetre excess 

There are large uncertainties in the physical properties of 
the VCGs. Such a large mass of dust locked up in very cold 
dust is not unlikely. However, there are some questions to 
address concerning the nature of these dust grains. How 
can a population of such very cold dust exist on galactic 
scales? How can this dust be so cold? Is very cold dust 
seen in other galaxies? In this section, we explore various 
hypotheses to explain the millimetre excess. 

We can also invoke a non-dust effect (i.e. molecular 
line and/or radio continuum contamination) or a cross- 
calibration problem. The non-dust effects have already 
been taken into account (see Sect. l2.1JI . Moreover the SED 
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of II Zw 40 (Galliano et al. 2003) shows the same trend in 
the millimetre regime, perhaps even suggesting that this 
excess might be a general feature of starbursting dwarf 
galaxies. 

4.3.1. Very cold dust in the literature 

A very cold component has been suggested for the Galaxy. 
For example, Reach et al. (1995) found a component with 
T = 5 — 7 K in some Galactic continuum spectra. More 
recently, Boulanger et al. (2002) published COBE and 
Archeops data for the Galaxy that implies a ~ 5 K dust 
component (with [3 = 2). This result has not yet been in- 
dependently confirmed. This component of very cold dust 
appears to follow the HI gas, suggesting that it is associ- 
ated with the cold neutral medium. Dupac et al. (2002) re- 
ported PRONAOS observations of the M 17 complex and 
found dust as cool as 10 K in some cold clumps. Popescu 
et al. (2002) observed a sample of 63 late-type galaxies in 
the Virgo cluster using ISOPHOT. Although their wave- 
length coverage is limited, observing only up to 170 /xm, 
they find a cold dust component is present in most cases, 
with temperatures as low as 10 K in irregular and blue 
compact dwarf galaxies. 

James et al. (2002) calculated dust masses for galaxies 
from the submillimctcr emission, assuming that the frac- 
tion of metals incorporated in the dust in galaxies is a 
universal constant. For NGC 1569 they determine a dust 
mass of ~ 1 x 10 5 M©, which is similar to the mass that 
we derive (1.6 - 3.4 x 10 5 M Q ). 

Lisenfeld et al. (2002) published a dust SED for 
NGC 1569 computed with the standard DBP90 model. 
They found no submillimctre excess and discount the pres- 
ence of very cold dust since they were able to explain the 
submm emission with the VSGs. We now discuss the rea- 
sons for our different results. We find that the use of the 
CVF spectrum in the MIR provides an important and 
tight constraint on the slope of the VSG size distribution 
in our model. Lisenfeld et al. (2002), however, characterise 
the MIR wavelength regime with only a single 12 \im IRAS 
data point. This effects the fit significantly. 

In addition to the different observational constraints 
used to construct the observed SED, Lisenfeld et al. (2002) 
used the DBP90 model in a very different, more limited 
way than that presented here. Since we have more data 
to constrain the DBP90 model, we allow the grain prop- 
erties to differ from those of the original DBP90 model, 
which was constructed to explain the Galactic dust emis- 
sion. Lisenfeld et al. (2002) use the original Galactic dust 
parameters to obtain their fit. Their input ISRF is that of 
the Galaxy, which is much softer than that of NGC 1569 
(see section ^21, and is scaled up by a factor of ~ 60. 
They did not investigate the dust size distribution param- 
eter space. Our results are based on a % 2 evaluation to ob- 
tain the best fit. Lisenfeld et al. (2002) conclude that the 
VSG component alone accounts for the FIR and submm 
and mm continuum. In our model, such a solution was 



attempted, but resulted in a very poor \ 2 since the MIR 
observations could not be well-fitted. 

4.3.2. Optical properties of BGs 

We have investigated other means to increase the emis- 
sivity in the submm, instead of invoking a VCG compo- 
nent. A change in the emissivity slope could be due to a 
change in the optical properties of the big grains. Indeed, 
Mennella et al. (1998) measure the absorption coefficient 
per unit mass (Q a bs) of cosmic dust analog grains, over 
the temperature range 24 — 295 K. The laboratory exper- 
iments, show a decrease of the emissivity index, /3, as the 
temperatures increases, effecting the Q a bs for wavelengths 
A > 30 /im. Our big grains reach higher temperatures 
than those of the Galaxy (Table HaJ) and, in the standard 
model, Q® b G s oc A -2 for A > 100 fj,m. We attempted to 
roughly reproduce the effect observed by Mennella et al. 
(1998) using QBg oc A~t for A > 100 /im with 7 < 2. 
However, this modification failed to give an acceptable fit 
to the data in the submm/mm wavelength range. 

Grain-grain coagulation, leading to porous aggregates, 
can increase dust emissivities at submm wavelengths. 
This has been invoked to explain the elevated opac- 
ity of circumstellar dusty envelopes and dense molecular 
cores (Ossenkopf & Henning 1994; Stepnik et al. 2003). 
However, this effect is localised and has not yet been seen 
on galactic scales, as would be necessary to explain the 
data presented here. 

4.3.3. Embedded dust as the source of the mm excess 

The properties of the dust differ considerably from one 
region to another in a galaxy, depending on the nature of 
the environment (HII regions, cirrus, dense core, etc.). The 
DBP90 model has been constructed to describe the diffuse 
ISM of the Galaxy. Consequently, applying it to an entire 
galaxy like NGC 1569 is not straightforward. However, 
NGC 1569 is a peculiar galaxy: it is a starbursting dwarf 
and its emission properties show that it resembles a gi- 
ant HII region (particularly in the MIR; see Madden et 
al. 2003). Moreover, the emission from hot dust is en- 
hanced relative to the cold dust in the starburst regions 
(see Fig. |2J. This departure could explain our inability 
to fit the SED with only the three standard components. 
The various phases of the ISM could have very different 
temperatures and technically should not be treated as if 
there were the same physical region. 

If we assume that the very cold dust is made up of 
grains deeply embedded in clouds, where the UV-optical 
radiation can not easily penetrate, the dust FIR-mm emis- 
sion could be a heating source for this dust component 
(Fig. I16fl . We can make a rough estimate of the temper- 
ature of the very cold grains (i.e, the VCGs) to verify 
this explanation. Let us assume that the photons which 
heat the VCGs are due to the emission from the BGs, 
simply modeled as a blackbody, made of grains of ra- 
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Fig. 16. Schematic diagram of a clump. 



dius a B Q, temperature T B g = 57 K (corresponding to 
the peak wavelength of the BG component), mass Mbg 
and emission coefficient Q om (aBGj ^bg), ^bg being the 
frequency where the maximum emission by the BGs oc- 
curs. These are the parameters that we derived from our 
fitting of the SED of NGC 1569. The VCGs are mod- 
eled with a blackbody of temperature Tvcg, mass Mycc 
and absorption and emission coefficients Q a bs(avCGi ^bg)> 
Qem(avCG, ^vcg), ^vcg being the frequency where the 
maximum emission by the VCGs occurs and are made 
of grains of radius aycc embedded in a clump of ra- 
dius R c . The luminosities emitted by these two compo- 
nents, Lbg an d LycG, can be determined from the num- 
ber of grains of each species, N x = 3M x /47ra x/ o x (x be- 
ing either BG or VCG and p x being the specific mass 
densities), and the luminosity radiated by only one grain 
C K = 47ra x Q om (a x , ^ x )ctT 4 . The total luminosity is then 
L x = N x £ x (x=BG or VCG). The ratio 6 = L B g/L V cg 
is 

/ M B g \ ( Pvcg\ ( avcc 
VMvcg/ 



e 



Pbg 

Hem (aBG,"BG) 



aBG 
Tbg 



(8) 



,Qem(avCG, ^VCG) / \ TyCG , 

and its numerical value from the model (Fig. [SJ) is O = 
4.1 x 10 4 . Using the same method, the heating and cooling 
luminosities for the VCGs are: 



J hcat 



-'cool 



3M B G 2 o i \ rrvi 

7ra BG y om ( i aBG, ^bg )® Ibg 



47ra| G p B G 
3M VC G 

47ra VCG^VCG 



47ra VCG Q em (avcG, *Aacg)o'T 



4 

VCG 



(9) 

which is equivalent to Lhoat = 1/4Lbg and L coo i = Ly GG . 
This factor of 1/4 is due to the fact that BGs do not 



radiate only in the direction of the VCGs inside the clump, 
they also radiate away without intercepting the VCGs. 
The powers absorbed and emitted by one VCG grain are: 



abs 



at 47ra VCG Q a b s (avcG, ^bg) 



4ttR2 

Lcool 



(10) 



4W 4 ™ VCG 



Equating absorption and emission gives: 

1 / Mbg \ / aycG \ / Pvcg \ 

-VCG 



4 V m vcg/ \ a B c / \ Pbg / 

Qcm(aBG, ^BG)Qabs(aycG, ^Bg) 
Qem(avCG, ^VCg) 



-,1/4 



T 



BG • 



(11) 

In the Eq.^2 the two unknowns are Tyco and aycG- We 
can deduce aycG by coupling Eas.lSlandlTTI 



6 ^ 4/Q abs (a V cG, fBc) 



(12) 



To be consistent with our model, which assumes (3 = 1, 
the VCGs are carbonaceous, pvcc = 2.2 gem -3 and 
Pbg = 3.3 gcm~ 3 . To determine the order of magni- 
tudes of the efficiencies, we refer to Draine & Lee (1984). 
For the BGs we use the astronomical silicate values at 
A = 90 /i 77i, Qabs/a — 0.017 fan . This is in the 
regime where Q a b s /a does not depend on the grain ra- 
dius. For the VCGs, we used graphite values. The coef- 
ficients are: Q a bs/a — 0.022 /im" 1 at A = 90 [rnn and 
Q abs /a ~ 4.5 x 10~ 4 [ivoT 1 at A = 660 fj,m. Eg. IT21 gives 
us the radius of the VCGs, aycc — 4 nm is also in the 
same regime where Q a bs/a is independent of a. With this 
simple analytical model, the temperature of the VCGs is 
estimated to be 6 K < Tyco < 9 K. Thus, the hypoth- 
esis of very cold dust hidden in clumps bathed by the 
FIR radiation from warmer grains is consistent and gives 
the appropriate order of magnitude of Tycc- Moreover, 
Bernard et al. (1992), using a radiative transfer model 
coupled with DBP90 in the Galactic case, showed that 
the temperature of the dust deeply embedded in a dense 
cloud is ~ 6 K. Radiative transfer may indeed turn out 
to be the key to understanding this property of the dust 
spatial distribution in these galaxies. 

We can attempt to estimate the clump radius (R c ). 
First, we can determine a minimum clump radius (RJ. nin ) 
by estimating the optical depth required to shield the stel- 
lar radiation, since in the above scenario, the VCGs are 
heated only by the envelope of BGs. We consider that 
the VCGs are screened by the BGs when the stellar en- 
ergy seen by the VCGs is at most equal to the energy 
emitted by the VCGs. Given our SED and the extinction 
curve of Fig. ^1 the optical depth in the V band needed 
is Ty hleld ~ 60. The average wavelength of the stellar ra- 
diation is Aisrf — 0.2 /im. At this wavelength, the optical 
depth needed is T shleld (i^isrf) — 200. At the wavelength 
where the BGs emit (A = 90 /im), the optical depth is 
t = 0.1, thus we verify that the shell is optically thin to 
this radiation. The contrast density between the clumps 
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(index c hereafter) and the intcrclump media (index icm 
hereafter) is a = n c /nj cm . The minimum radius of a clump 
R™ m is the radius required to reach the optical depth of 
T shieid^ v ^^y If we assume that only the BGs are respon- 
sible for the shielding, 

T shield (vi SK F) = RfV G ™i G Q abs (a BG , ^isrf) (13) 

where n^ G is the density of the BGs in the clumps and 
Qabs(aBG> ^isrf) — 0.1 from Draine & Lee (1984). The 
actual clump size is difficult to estimate without some as- 
sumptions for n|? G or a. Let us assume that n|? G is the 
homogeneous density n^ G n ~ 3 x 10~ 8 cm~ 3 scaled by 
a. From Eq. H3 we find that R™ in ~ 8 pc if a = 10 4 , 
the extreme case, where there is only one clump, and 
R™ in ~ 0.8 pc if a = 10 5 . To estimate the size of the 
VCG core, we require that it is optically thin to the BG 
radiation. The ratio between the mass of BGs and VCGs 
constrains the density of this core. The maximum ra- 
dius of the VCG core is ~ 4 — 6 pc when a = 10 4 , the 
number of clumps is 1 or 2 and the volume filling factor 
4> ~ 3 — 6 x 1CP 5 . When a = 10 5 , the radius of the core is 
~ 0.4 — 0.6 pc, the number of clumps is between 200 and 
300 and the filling factor ~ 3-6 x 10~ 6 . This rough esti- 
mation tends to imply that there would be a small number 
of dense clumps in this galaxy. Berkuijsen (1999) gives a 
range of 10~ 3 — 10~ 5 for the H2 volume filling factor of 
the Galaxy which is greater or equal to the range we find 
in NGC 1569. This clumpy characteristic is also consis- 
tent with relatively bright 158 jim [CII] emission seen in 
NGC 1569 and in low-metallicity environments, in general 
(Jones et al. 1997). 

5. Summary and conclusions 

Our modeling of the dust SED in NGC 1569 suggests that 
the nature of the dust in low-metallicity environments dif- 
fers from that of metal-rich galaxies. 

1. We have presented new SCUBA images at 450 and 
850 fim. With additional data from the literature we 
constructed the observed dust SED for this galaxy. 

2. We have used a stellar evolution model (PEGASE), a 
photoionisation model (CLOUDY) and a dust model 
(DBP90) to compute a self-consistent theoretical dust 
SED, independently synthesizing a global ISRF from 
UV-to-optical observations with further constraints 
from MIR- ionic lines. The DBP90 model has been used 
to investigate a large range of dust parameters, mak- 
ing the link with the extinction deduced from the dust 
model and the deredenning of the UV-to-optical data. 

3. The synthesized global ISRF is consistent with a very 
young population produced by the recent starburst 
and an older population tracing the underlying stars 
of the galaxy. 

4. We find very low abundances of PAHs, and smaller 
overall sizes of grains emitting in the MIR and FIR, 
compared to the Galaxy. Due to the small sizes, most 
of the grains are in a stochastic heating mode, and not 



in thermal equilibrium with the radiation field - even 
grains emitting at FIR wavelengths. 

5. The presence of a millimetre emission excess can be 
explained by ubiquitous clumps of very cold (5 to 7 K) 
dust grains with (3=1. These very cold grains can be 
responsible for 40 to 70 % of the total dust mass in the 
galaxy, M dust = (1.6 - 3.4) x 10 5 M Q . The absence of 
important submm/mm observations can result in, at 
least, an order of magnitude of dust mass being missed. 

6. The gas-to-dust mass ratio ranges from 740 to 1600, 
greater than that in the Galaxy, even taking into ac- 
count the relatively large mass of very cold grains. 
The dust-to- metals mass ratio ranges from 1/4 to 1/7, 
smaller than the Galactic value. 

7. From the dust size distribution, we generate an extinc- 
tion curve for NGC 1569 (which has similar character- 
istics to that of the LMC). We also derive a synthesized 
radiation field and a wide range SED (UV - mm), for 
the global galaxy. 

8. The bulk of the cold dust traced by the submm/mm 
observations is distributed between star forming re- 
gions, not concentrated toward the outer regions of 
the galaxy. Our results are consistent with a clumpy 
medium and a filling factor lower than 10~ 4 . 

Due to the low metallicity of NGC 1569, the ISM would 
appear to be very clumpy, and this has important conse- 
quences for the dust models. From our detailed dust mod- 
eling of dwarf galaxies (see also Galliano et al. 2003), we 
caution that dust models assuming the Galactic dust size 
distribution, Galactic extinction curve and Galactic ISRF 
will be in error. Models not dealing with stochastic heating 
processes should be dealt with cautiously, as this is prob- 
ably the most important physical process for dust heating 
in dwarf galaxies. Another caution that follows from this 
work is the importance of obtaining data longward of the 
FIR. This is presently easier said than done. However, in 
the near future we will have the capability to sample the 
submm/mm emission from galaxies in more detail, and 
with greater sensitivity, with SOFIA, ASTRO-F, Herschel, 
Planck and ALMA. Without a proper understanding of 
the submm/mm SED of galaxies, galaxy number counts 
could be incorrect since the excess submm/mm emission 
might erroneously be attributed to higher star formation 
activity in redshifted galaxies, instead of the presence of 
a large, very cold dust mass. 
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